Mon. Not. R. Astron. Soc. 000,[T]-?? (2013) Printed 18 January 2013 (MN MgC style file v2.2) 



The host-galaxy response to the afterglow of GRB 100901A* 



Olga E. Hartoog^'^ t, Klaas Wiersema^, Paul M. Vreeswyk"^, Lex Kaper\ 
Nial R. Tanvir^, Sandra Savaglio^, Edo Berger^, Ryan Chornock^, Stefano Covino^, 
Valerio D'Elia^'^, Hector Flores^^, Johan R U. Fynbo^\ Paolo Goldoni^^, An- 
- - - dreja Gomboc^^'^'^, Andrea Melandri^'^^, Alexei Pozanenko^^, Joop Schaye^, Anto- 
^ nio de Ugarte Postigo^^, Ralph A. M. J. Wijers^ 



O 



O 

u 

o 
is 



> 

(N 

rn 

O 

m 



X 



(Affiliations can be found after the references) 



18 January 2013 



ABSTRACT 

For Gamma-Ray Burst 100901 A, we have obtained Gemini-North and Very Large Tele- 
scope optical afterglow spectra at four epochs: one hour, one day, three days and one week 
after the burst, thanks to the afterglow remaining unusually bright at late times. Apart from 
a wealth of metal resonance lines, we also detect lines arising from fine-structure levels of 
the ground state of Fe ii, and from metastable levels of Fe ii and Ni ii at the host redshift 
(z - 1.4084). These lines are found to vary significantly in time. The combination of the data 
and modelling results shows that we detect the fall of the Ni ii"*/^9/2 metastable level popula- 
tion, which to date has not been observed. Assuming that the population of the excited states 
is due to the UV-radiation of the afterglow, we estimate an absorber distance of a few hundred 
pc. This appears to be a typical value when compared to similar studies. We detect two inter- 
vening absorbers {z - 1.3147, 1.3179). Despite the wide temporal range of the data, we do 
not see significant variation in the absorption lines of these two intervening systems. 

Key words: gamma-ray burst: individual: GRB 100901A; gamma-rays: bursts - galaxies: 
abundances - galaxies: ISM - galaxies: distances and redshifts 



1 INTRODUCTION 

Shortly after the first detection of an optica l afterglow asso- 
ciated with a long gamma-ray burst (GRB, Ivan Paradiis et alj 
Il997h . it became clear that GRBs may be useful probes of 
distant galaxies: afterglows can be very bright and have a simple 
power law continuum at ultraviolet (UV) to optical wavelengths, 
against which otherwise undetectable absorbing systems, such 



(see 



Vreeswiik et al 



Vreeswii k et al. 2001; 

Berger et al. 2006; 



Prochaska et al 



12004; 



Fynbo et al 





2003 




2006 



l2007h . Long G RBs are formed as end products 

of the evolution of massive stars dWooslevIl 1 9931 ; IPaczvnskill 1 9981) 
and should therefore probe star-forming regions. Their optical 
afterglows fade away on short timescales (~days), making it 
possible to study their host galaxies in emission as well (see e.g., 



* Based on observations obtained at the Gemini Observatory at Mauna Kea 
under programme GN-2010B-Q-7, and observations obtained with ESO 
Telescopes at the Paranal Observatory under programme 085.A-0009(B). 
t 0.E.Hartoog@uva.nl 



ISavaglio. Glazebrook & Le Borgnell2009l ; iHiorth et al.ll2012h. For 
a review of long GRBs and their host galaxies, see e.g., Savaglid 
l l2006Ll2012h . 

Optical spectra (i.e., rest frame UV at redshifts around z ~ 
1 - 2) of GRB afterglows generally show strong absorption lines 
of metals, both at low and high-ionisation stages, as well as 
high neutral hydrogen column densities (see Ipvnbo et al.l l2009l 
for a large and relatively unbiased sample). Absorption lines pro- 
vide the redshift and allow estimation of column densities of the 
host galaxy of the GRB and possibly other absorbers along the 
line of sight. This provides valuable information on the metal- 
licity and dust content of galaxies at arbitrary redshift, regard- 
less of the brightness of these systems, which is very difficult 
to obtain otherwise. The highest redshift lines in the afterglow 
spectrum formally only provide a lower limit on the host red- 
shift. However, the presence of fine-structure lines can confirm 
the galaxy as the host, because these arise from excited states that 
are expected to have been populated by the UV-afterglow radia- 
tion o f the GRB jProchaska. Chen & Bloonj200d : IVreeswiik et all 
|2007|) . These lines are not seen in foreground line-of-sight ab- 
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sorbers in quasar (QSO) spectra, except for some relatively low- 
energy excited states from carbon. Due to the transient nature 
of the afterglow, the absorption lines from these excited states 
are expected to vary in strength, which may provide informa- 
tion on the internal distribution of gas and kinematics within the 
host galaxy; often the only way to obtain such information for 
these distant and mostly faint galaxies. Fine-structure line vari- 
ation has been measured and modelled for a handful of bursts 



Ce.s.. Vreeswiik et alj 20071: Dessauges-Zavadskv et alj 20061: 


D'EUa et al.l 2009al: iLedoux et al. 


l2009l: 


de Uaarte Postieo et alj 


201 ll: De Cia et all |2012|: Vreeswiik et al. 


2013f). which allowed 



constraints to be placed on the distance betwe en the burst and the 
absorbing material (see IVreeswiik et al.ll2013l for a description of 
the methods used). 



The primary difficulty in probing column density variations in 
GRB sightlines to fit to the models described above, is the require- 
ment on afterglow brightness: for fainter afterglows the necessary 
high signal-to-noise and high spectral resolution are difficult to ob- 
tain with current instrumentation. This also limits this method to 
early-time data and a limited range of instruments, making it diffi- 
cult to build up a sample. In the cases described above, the variation 
is measured with high resolution spectrographs, and all spectra are 
obtained within a few hours post burst (observer's frame). 

The optical afterglow of GRB 100901 A remained unusually 
bright out to very late times, which allowed us to collect spectra 
over a time span of a week with the Gemini-North Multi-Object 
Spectrograph and X-shooter on the ESO Very Large Telescope. In 
ffiis paper we investigate the behaviour of ffie excited level popu- 
lations over this unusually long time span, and the application of 
the excitation models to low and intermediate resolution spectra. 
The observations and data reduction are described in Section [2] 
We study the metal resonance lines at the host-galaxy redshift in 
Section [3711 In Section [3. 1.1 1 we examine the dust content of the 
host galaxy. The two intervening absorbers are analysed in Sec- 
tion [3T3] In Section [3j2] the fine-structure line variability is de- 
scribed and modelled. In Section [331 we briefly discuss the emis- 
sion lines from the host galaxy. We discuss our results in Section|4] 
and summarise our conclusions in Section|5] 

Throughout the paper we adopt a standard ACDM cosmology 
with Ho = 71 kms"' Mpc"', fi„ = 0.27, Qa = 0.73. We indicate 
line transitions by their vacuum wavelengths. 



2 OBSERVATIONS 

GRB 100901Awas detected by Swift on 2010 Se ptember 1 at 
13:34 : 10 UT. With a total dur ation of Tga = 439±33 s jlmmler etal] 
l20ld : ISakamotoet aDl20IOh the burst is clearly classified as a 
long burst. The optical afterglow candidate at RA (J2000) = 
0lH9"'03.42% Dec =+22°45'30.8" (90% confidence error radius 
of about 0.81") was identified in an Ultraviolet/Optical Telescope 
(UVOT) exposure w hich started 147 s a fter the Burst Alert Tele- 
scope (BAT) trigger dlnnmler et al.ll2010l) . Automatic observations 
with the Faulkes Telescope North identified an unca talogued object 
at the position consistent with the UVOT candidate jGuidorzi et alj 
I2OIOI) . Several follow-up photometry efforts followed, resulting in 
the light curve presented and analysed in Gomboc et al. (in prep) 
and iGorbovskov etal. ( 20121) (see also Figure [8). The first opti- 
cal spectrum of the afterglow was taken with the Gemini Multi- 
Object Spectrograph (GMOS) on Gemini-North (programme GN- 
2010B-Q-7, PI Ta nvir), approximately I hr 15 mins (0.0526 days) 
after burst trigger (IChomock et al .l2010b . With this instrument, two 



other spectra at respectively 1 and 7 days after the burst were ob- 
tained. A fourth spectrum was obtained approximately 3 days af- 
ter the burst with the X-shooter spectrograph mounted on the ESO 
Very Large Telescope (VLT), under programme 085.A-0009(B) (PI 
Fynbo). See Table[T]for a detailed log of the observations presented 
in this paper. 

2.1 Gemini-N/GMOS spectroscopy 

Gemini-N/GMOS is a low-resolution long slit spectrograph, 
equipped with three detectors. For the observations presented here, 
a slit width of 0.75" and the B600 grism with the G5307 order 
suppression filter have been used, which resulted in the wavelength 
coverage and resolving power reported in Table [T] We used four 
exposures per epoch, using dithers in both dispersion and spatial 
coordinates to sample over the chip gaps and regions affected by 
amplifier location. The GMOS spectra have been reduced with the 
IRAF packages for Gemini GMOS (version 1.9), using arc line ex- 
posures taken directly before and after the science data. The four 
exposures are combined after extraction. The resulting spectra were 
normalised and no flux calibration was performed. 

2.2 VLT/X-shooter spectroscopy 

VLT/X-shooter is a cross-dispersed echelle spectrograph where the 
incoming light is split into three wavelength arms using dichroics, 
covering the full optical to near infrared wavelengt h range simulta- 
neously to'Odorico et al.ll2006l : IVemet et al]|201 ih . The X-shooter 
spectra have been taken in nodding mode with 1x2 binning (i.e., 
binning in the dispersion direction) in the UV-Blue (UVB) and Vi- 
sual (VIS) arms, using a lOOkHz / high-gain readout and 4 x 600s 
exposure times. The selected slit widths were I.O, 0.9 and 0.9" for 
the UVB, VIS and near-infrared (NIR) arms, respectively. A 5" nod 
throw along the slit was used to improve sky subtraction. The spec- 
tra were reduced using the ESO pipeline so ftware version 1.2.2 , 
using the so-called p hysical model mode i lGoldoni et al] l2006l : 
iModigliani et alj|2010h . We used calibration data (bias, dark, arc 
lamp, flat-field and flexure control frames) taken the next day or as 
close in time as possible to the science observations. The bin sizes 
of the resulting spectra were 0.2, 0.2, and 0.5Afor UVB, VIS and 
NIR spectra, respectively. Extracted spectra were flux calibrated 
using spectrophotometric standard star exposures, also taken the 
same night, resulting in roughly flux-calibrated spectra spanning 
the near-ultraviolet to the near-infrared. We caution that the weather 
conditions during this night were poor, and likely non-photometric. 
We used a telluric standard star (HD 4670, a B9 V star) to correct 
the telluric absorption feature s in the NIR afterglow spectrum, us- 
ing the IDL SPEXTOOL package jVacca. Gushing & Ravnerl2003l , see 
also lwiersemall201]] for the use on X-shooter data). 



3 ANALYSIS 

The GMOS epochi spectrum (see Figure[T) is that of a bright after- 
glow including many metal absorption lines. At a redshift of 1 .4084 
we detect strong resonance lines of Fe 11, Mn 11, Cr 11, Al 11, Al m, 
Znn and Sin. Additionally, at the same redshift, line transitions 
arising from excited states of Fe n and Ni n are clearly detected 
(see Sec tionl3.2[(. Therefore we identify this as the host-galaxy red- 
shift Zh jChomock et alj|2010l) . We do no t detect lines from F e m, 
as observed in the case of GRB 080310 jPe Cia et al.]|20I2l) . No 
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Figure 1. Normalised GMOS spectrum (epochl) of the afterglow of GRB 100901 A. The red labels indicate the absorption lines of the host {zi, = 1.4084); the 
green and blue labels (with vertical olfset) indicate the absorption lines formed in two intervening systems at zi = 1.3147 and Z2 = 1-3179. The magenta labels 
show unshifted Na i lines (z = 0), resulting from the absoiption in the Milky Way. The numbers at the labels indicate the configuration of the lower level and 
the vacuum rest wavelength of the transition. The host shows resonance lines from the ions Mn ii, Fe ii. Cm, Zn ii, Sin, Al ii and Al iii, as well as lines from 
excited states of Feii and Niii, which are found to vary significantly (see Section lX2t . In the intervening systems we detect Mgi, Mgii, Pen, Alii and Aim. 
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Table 1. Overview of the spectroscopic observations of the afterglow of GRB 100901A, arranged by instrument and epoch. Columns (1) and (2) indicate 
telescope, instrument, setup, and slit dimensions. Column (3) gives the wavelength coverage for this setup and the resolving power R = A/AA, with AA the 
FWHM of an unresolved line. Column (4) gives the mid time of the epochs of observation in days after the burst BAT trigger. Column (5) lists the interpolated 
observed R-magnitude of the afterglow at this epoch (Gomboc et al., in prep). Column (6) lists the total exposure time of the observation in this epoch. Columns 
(7) and (8) give the seeing and the aimiass. The seeing is measured by fitting a gaussian along the spatial direction of the slit, around the central wavelength of 
each spectrum. Column (9) gives the signal-to-noise ratio (S/N) in the continuum at 4150 and 6100 Aper GMOS pixel (6A = 0.91 A). Column (10) gives the 
name we will use to refer to the epoch. 



Telescope 


Instrument 


Spectral Range 


Time after Burst 


Mr 


Exp. Time 


Seeing 


Airmass 


S/N/(5.1 = 0.91A) 


Abbrev. 




and Setup 


and Resolution 


(days) 




(s) 


(") 




at 4150, 6100 A 




Gemini-North 


GMOS, 


3810 -6710 A 


0.0526 


17.8 


4x400 


0.62 


1.06 


38,95 


epoch 1 




B600+ G5307 


R= 1000- 1900 


1.0095 


18.6 


4x500 


0.52 


1.01 


16, 70 


epoch2 




0.75" X 330" 




7.0270 


22.3 


4x 1200 


0.68 


1.05 


5, 13 


epoch3 


Vei7 Large 


X-Shooter 


3000 - 25000 A 


2.7492 


20.2 


4x600 


1.30 


1.48 


9,7 


xsh 


Telescope 


UVB 1.0" X 11" 


R = 5400 


















VIS 0.9" X 11" 


R = 7400 


















NIR 0.9" X 11" 


R = 5800 

















Lyman-tt is detected, because it is too far in the blue for the cov- 
ered wavelength range. Two intervening absorbers are visible at 
zi = 1.3147 and zi = 1.3179 dChomock et"ai]|201(ll) . In these sys- 
tems we detect Fe n, Al n, Al m, Mg i and Mg ii (rest frame equiv- 
alent width W,\ < 1 A). We note that these lines are weaker at the 
intervening system redshifts than at the host redshift. There are no 
fine-structure lines detected at any of the intervening absorber red- 
shifts. Nai 5892, 5898 are detected at z = and are thus due to 
foreground absorption in the Galaxy. 

The GMOS epoch2 spectrum (not fully shown in a figure) 
looks similar, though the signal-to-noise ratio (S/N) is lower due 
to the decreasing brightness of the afterglow. The line transitions 
arising from the metastable Niii ^Fcj/t level (see Table lA4T l appear 
more pronounced in epoch2 than in epoch 1; the Fen fine-structure 
lines have become slightly weaker 

Since GMOS epoch3 is observed after the break in the light 
curve at ~ 2 days (Gomboc et al., in prep), which causes the bright- 
ness to drop more quickly, the S/N is lower than for the other 
epochs, despite the longer integration time and comparable weather 
conditions. Still, strong metal resonance lines are clearly detected. 
Fine-structure lines that were present in epoch 1 and lines from 
metastable levels that were present in epoch2 are much weaker or 
not detected any more (see Section|3j2}- 

The quality of the X-shooter spectrum (xsh, after 2.75 days, 
see Table [TJ is low due to poor weather conditions and large air- 
mass, but because of the extraordinarily large wavelength cover- 
age of this instrument, this spectrum does provide additional in- 
formation. Moreover, the higher resolution allow us to exclude the 
existence of different velocity components down to ~ 70 km s"' 
with respect to ~ 200 km s"' for GMOS. The strong and saturated 
Mg II .1.1 2796, 2803, Mg 1.12853, and Caii.i.i 3934, 3969 absorp- 
tion lines from the host galaxy are outside the range of GMOS 
due to the redshift, but they are clearly detected in xsh- VIS (not 
shown in a figure). Na 1.15892, 5897 falls in a region with atmo- 
spheric absorption and can therefore not be detected. We do not 
detect the red wing of the Ly-o- absorption line due to low S/N. 
Furthermore, in xsh-NIR, we marginally detect forbidden oxygen 
emission line [Om] A5007 at the host-galaxy redshift (see Figure |9] 
and Section [33] ). 

We do not detect prominent interst ellar (dust) extinction 
features such as the 2175A feature (see e.g., lElfasd6ttir et alj2009l) 
or dilfuse interstellar bands in any of the spectra at the host-galaxy 
redshift. 



3.1 Resonance lines 

To measure the observed equivalent width VKi.obs of the absorption 
lines in the spectrum, we locally fit gaussian functions to the nor- 
malised spectrum with ngaussfit, a task in the stsdas package in 
IRAF. In principle, Wj.obs can be obtained directly by integration, 
but as can be seen in Figure [T] many lines are blended. In some 
cases one line needs to be fixed in order to measure the potential 
variation of the other, ngaussfit can de-blend up to three lines, with 
the possibility to keep part of the parameters fixed. This method 
allows that in spectra with low S/N, an absent or very weak line is 
sometimes formally best fitted with a line with negative equivalent 
width (i.e., an emission line, see also Figure|6]and|7]l; for these lines 
we obtain upper limits for IV.i.obs. 

The formal error on the equivalent width is obtained from the 
error spectrum: 



crW, 



1 /■ 



(1) 



with 6A the spectral bin width, cr, the value of the error spectrum 
in bin ; and F,. the continuum flux. The summation runs over two 
times the fuU-width-at-half-maximum (FWHM). Furthermore, we 
add 3% of the equivalent width to the formal error to account for 
the uncertainty in placing the continuum. In the low-resolution 
GMOS spectra, all absorption lines are unresolved and can be 
fit with a gaussian with a fixed FWHM of 3.6 A (the resolution 
element), which is determined from the width of the lines in the 
arc-lamp frames. In xsh we fix FWHM = 1.0 A and 0.7 A in 
the UVB and VIS-arms, respectively, based on the widths of arc 
lines and telluric emission lines. Some strongly saturated lines 
(for example the Mg ii doublet) in these spectra cannot be fit with 
gaussians and are therefore integrated numerically. The resulting 
rest frame equivalent widths = W,i.obs/(l + a) of the resonance 
lines in the first two epochs are listed in Table |2] 

To infer column densities of the dilferent ions from the absorp- 
tion line equivalent widths, we make use of a m ulti-ion single com- 
ponent curve-of-growth (MISC-COG) analysis ( ISpitzeJl978l) . The 
Doppler line width is in principle related to the temperature and the 
level of turbulence in the absorbing gas. Limited by spectral reso- 
lution, it is difficult to distinguish between this and broadening due 
to the complex physical structure of absorbing clouds in the sight 
line. Therefore we speak of an 'effective' Doppler parameter, which 
does not carry information about the temperature of the gas, but is 
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merely a simplification in the model jjenkinslflgseh . The assump- 
tion we need to make in this analysis is that the effective Doppler 
parameter b is the same for all ion species, and that it is single (one 
velocity component). In other words, we assume that the velocity 
structure is the same in all lines, and that it is dominated by small- 
scale effect such as thermal velocities and turbulence and not by 
differences in bulk velocities of different gas clouds. 



3.1.1 Host galaxy column densities 

We use three unblended Feii resonance lines (/12249, /12260 and 
/12586) to determine b for the line-of-sight gas in the host galaxy 
(see Table O. Because these lines cover both the linear and the 
flat part of the COG, both b and the column density for Fe ii, A'peii, 
can be constrained. They are dependent variables, therefore fit them 
simultaneously with a grid method, yielding b = 22.1+j g km s"' 
and logA'peii/cm-^ = 15.23!°;^^. 

The COG with the best fit value for h is shown in Figure |2] 
The ion species other than iron have line transitions either on the 
linear or on the flat part of the COG; therefore, b has to be assumed 
to be the same for all ion species in order to determine the column 
densities. The resulting column densities of the ions are listed in 
Table [2] For most lines we have taken an average of the first two 
epochs, except Mgii and Can which can only be measured in 
xsh. The equivalent widths of the Al 11 lines are slightly larger in 
epoch 1 while those of Aim are larger in epoch2. While not very 



significant, it may b e a sign of io nisation due to the GRB afterglow 
(see e.g., Vreeswii k et alj|2013h . Due to the low significance we 
decided not to take this into account in the excitation modelling in 
Sectionim 



Spectral resolution limits the ability to distinguish different ve- 
locity components in absorption lines: an instrumental FWHM = 
3.6 A (CMOS spectra) corresponds to ~ 200 km s~'. In xsh- VIS, 
a velocity component difference of 70 km s"' would be detectable, 
but we do not see indications for more than one velocity component 
in any of the lines. This could be due to the low S/N of xsh, because 
most intermediate to high-resolution spectra of GRBs show absorp- 
tion l ines with at least two components (but see e.g.. lLedoux et alj 
I2OO9I) . 

IProchaskal (l2006h warns that column densities derived with 
the MISC-COG method are systematically underestimated, espe- 
cially when the underlying line structure is that of separated clumps 
of gas at different velocity shifts. In this case one derives a high b 
with MISC-COG analysis, and therefore a lower column density. 
Results for a single component C OG analysis with b > 20 km s"' 
are highly suspect, according to IProchaskal (l2006h . That COG 
analysis of a low-re solution spectrum can be accurate is shown by 
iD'Elia etal.l l l20Ilh . who compare low and high-resolution spectra 
for the afterglow of GRB 081008. The COG analysis gives the 
same result as line fits to the high-resolution spectrum within 3cr. 
They conclude that this is linked to the low level of saturation 
of the lines. The equivalent widths of the low-ionisation lines in 
GRB 100901A are found to be in the lo west 10% of the sample 
of 69 low-resolution spectra conducted by Ide Ugarte Postigo et alj 
( l2012l) . which strengthens the case that a COG analysis on this 
spectrum can provide accurate results. However, column densities 
resulting from COG analysis always need to be treated with 
caution. 




Figure 2. Curve-of-growth (COG) for the resonance lines in the host 
galaxy of GRB 100901A. The solid line gives the best fit COG with b = 
22. 1 km s"' , the dotted lines show the effect of the en'or on b. The dashed 
line, which coincides with the lineal' part of the COG, is a COG for = 00, 
that can be used as an approximation for weak lines. 
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Figure 3. [Fe/Zn] as a function of the Znn column density in QSO- 
DLAs (red diamonds), GRB-DLAs ( black circles) and GRB 100901A(blue 
square). QSO-DL A values a re from Luetal. (1996); Prochaska & Wolfg 
(l996. 1997); Pe ttini et all (_1999, 2000); Prochaska & Wolfe (200(|^ 



[Prochaska et alj i2 0()l). GRB- DLA values are fromlSavaglio. Fall & Fiord 
)2003l): iFvnbo etal (2002): Ipiore et al.l ll2005h: IWatson et al. ( 2QQ(| 
IPenprase et al.l <2006l) : lvreeswiik et al.l ( |2007|) : Id Elia et al.l <2009bll201 11) 



3.1.2 Dust depletion 

The strength of the absorption lines in the spectnam are a measure 
of the abundance of an element in the gas phase only. Apart from 
the overall gas-to-dust ratio, the fraction of a specific element that 
is locked onto dust depends on the element species and on the na- 
ture of the dust. Therefore, the pattern of relative abundances, the 
depletion pattern, might give insight into the amount and nature of 
the dust. In Table[2]we list the abundance ratios with respect to the 
refractory (depleted) element Fe and the non-refractory element Zn, 
compared to the ratio in a solar abundance environment. We use the 
following, commonly used notation 
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Table 2. The measured rest-frame equivalent width W,i for the resonance lines at the host redshift, and the corresponding ion column densities, obtained 
with curve-of-growth (COG) analysis. The first column lists the ion and the vacuum rest wavelength (A); the second column gives the oscillator strength /j. 
Columns (3) and (4) give the rest-frame W,i in A measured in epochl and epoch2 (GMOS) respectively. For lines for which we have a good measurement both 
in epochl and epoch2, the weighted average is used; otherwise we use the value measured in epochl. The values that are finally used for the COG analysis 
are listed in Column (5). Column (6) gives the best fitting column density for the ion, fitted to a CO G with effective Dopp ler parameter b = 22.1^| g km s"' . 
Columns (7) and (8) give the relative abundances with respect to iron and zinc using solar values from lAsplund et all ( |2009|) (for notation see Equation|2)- The 
references for the wavelength and oscillator strength are in column (9). 



Line /, W,i/A Wa/A < WJA > \og(N-^^Jcm-^) [X/Fe] [X/Zn] ref 

epochl epoch2 

Aln 1670.7874 1.740 0.418 ± 0.093 0.338 ± 0.098 0.380 ± 0.067 13.57+° j\ -0.61 ±0.42 - 1.84 ±0.42 1 

Al in 1854.7164 0.539 0.113 ±0.054 0.162 ± 0.060 0.135 ±0.040 U.02"^° \l " " 1 

Aim 1862.7895 0.268 0.044 ±0.051 0.113 ± 0.056 0.075 ± 0.038 ... 1 

Sin 1808.013 2.080x 10-' 0.292 ± 0.062 0.267 ± 0.068 0.281 ± 0.046 15.96+g}^ +0.72±0.19 -0.51±0.19 2 

Znn 2026.136* 5.010 xlO"' 0.362 ± 0.046 0.280 ± 0.053 0.327 ± 0.035 13.52+555J^ +1.23 ±0.10 3 

Znn 2062.664 2.460x 10"' 0.191 ±0.041 0.206 ± 0.048 0.197 ±0.031 ... 3 

Crn 2056.2539 1.030x 10"' 0.150±0.040 0.151 ±0.046 0.151 ± 0.030 13.73+55°^ +0.36±0.12 -0.87±0.12 3 

Crn 2066.161 5.120x 10"^ 0.096 ± 0.038 0.111 ±0.045 0.102±0.029 ... 3 



Fen 2249.8754^ 


2.190 X 


10-3 


0.129 ±0.034 


0.151 ± 


0.040 


0.139 ±0.026 


15, 


91 +0.07 
■^•'-0.08 




- 1.23 ±0.10 


4 


Fen 2260.7793 


2.620 X 


10-3 


0.152 ±0.034 


0.168 ± 


0.040 


0.159 ±0.026 








4 


Fen 2586.6495 


7.094 X 


10-2 


0.705 ± 0.044 


0.727 ± 


0.048 


0.715 ± 0.032 










4 


Mnn 2576.877 


3.610 X 


10-' 


0.265 ± 0.032 


0.257 ± 


0.037 


0.262 ± 0.024 


13, 


90+0.04 
•^-0.05 


+ 0.07 ± 0.09 


- 1.16 ±0.04 


5 


Mnn 2594.499 


2.800 X 


10-' 


0.206 ± 0.030 


0.207 ± 


0.034 


0.207 ± 0.023 








5 


Mnn 2606.462 


1.980 X 


10-' 


0.164 ± 0.064 


d 




0.164 ±0.064 










5 


Mgn 2796.352 


6.123 X 


10-1 








0.920 ± 0.147'" 


15, 


Qg+0.83 
■""-0.60 


- 0.27 ± 0.83 


- 1.50 ±0.83 


6 


Mgn 2803.531 


3.054 X 


10-1 








0.918 ±0.126" 








6 


Can 3934.777 


6.500 X 


10-' 








0.907 ± 0.220'' 


13, 


74+0.47 
■'^-0.50 


- 0.33 ± 0.50 


- 1.56 ±0.50 


7 


Can 3969.591 


3.220 X 


10-' 








0.836 ±0.185'' 








7 



"The quantities [Al in/Fe] and [Al in/Zn] do not have any physical meaning, since the Al iii is probably mainly from a different region than Al n, Fe n and 
Znii (see also^avagli^^^all^004). 

*The contamination by Cr n /12026 is determined from the other Cr n lines and is neghgible. The contribution of Mg i ^2026 cannot be determined, but we 
assume it is small based on the absence of Mgi/11827. 

'Fen 2344, 2374, 3282 and 2600 are also clearly detected, but those are blended with fine-structure lines and therefore not used in the abundance study. 
Estimates of their equivalent widths can be found in the appendix. 

''The equivalent width of the line is fixed in the second epoch in order to measure the varied strength of the blended fine-structure line(s). 
■^Absorption lines are outside the GMOS range, measurement is from xsh spectrum. 

References for atomic data; (1) NIST Atomic Spectra Database, (2)^ergeson^Lawler Q993b), (3)^e2eson^Lawler (T993a) (4)^meretal^Q999^, 
(5) Morton (700J), (6) Vemer et al. (1996') . (71 Morton (19911 . 



[X/Y].logf™l-logfiil^) , 



with A' (X) the column density of element X, and 1 1 (X) the num- 
ber d ensity of element X in a solar environment dAsplund et al.l 
I2OO9I) . Differences from zero indicate dust content. We thus as- 
sume that relative abundances (both in gas and dust) are solar. Fe, 
Zn and Cr are iron peak elements, which makes this a reasonable 
assumption, alt hough the main product ion pathway of Zn is not 
fully identified dUmeda & Nomotol2002h . Discrepancies in relative 
abundances might also be due to a different star-formation history 
than the Milky Way. 

In Figure[3]we show [Fe/Zn] as a function of the Zn column 
density for the host of GRB 100901A, compared with a sample 



of other GRB-DLAs and QSO-DLAfl The relative abundance 
values indicate a stronger depletion (i.e., more dust) in GRB-DLAs 
than in QSO-DLAs, but among the GRB-DLAs, GRB 100901 A is 
not a particularly special case. Both GRB-DLAs and QSO-DLAs 
appear to follow the overall trend that systems with stronger Zn 
column densities show stronger depletion. The fact that the two 
classes populate different regions in the diagram is mainly a sight 
line effect. QSOs probe galaxies in random orientations with only 
a small chance that the densest region of a galaxy is intersected, 
while sight lines to GRB afterglows probe the gas towards the pro- 
genitor r egion, so they will in general show high column densitie s 
(see also lSavaglidl2006l ; |Prochaska et al.ll2007l;lFvnbo et al.ll2009h . 



DLA: Damped Lyman-o- system; a sight line absorber with A'hi > 2 x 
10^° cm-2 ^olfe. Gawiser & Prochaska 2()o3). A DLA system in a GRB 
afterglow spectrum is usually due to the host galaxy. 
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Furthermore, for GRB-DLAs with relatively low column densities 
and weak absorption lines, [Fe/Zn] is not well constrained (i.e., 
not reported in literature) because the fading afterglow prevents 
building up S/N over several nights of observations. This is not the 
case for QSO-DLAs, which is why these objects also populate the 
low A^znii part of the parameter space. 

A more detailed analysis of the dust content can be per- 
formed by comparing the heavy-element dust depletion pattern of 
GRB 100901 A wit h observations in the inters tellar medium (ISM) 
of the Milky Way ^Savage & Sembachl ll993) . We use the method 
described inlSavaglio & Fall (2004). We consider the observed de- 
pletion patterns in the MW sight lines trough a warm halo (WH), 
warm disk and halo (WDH), warm disk (WD) and cool disk (CD) 
as models. For the fit, we use two free parameters. One is the dust- 
to-metal ratio relative to the Galactic values k = koj^^/kj (J is one 
of the 4 depletion patterns). The second is proportional to the metal- 
licity (not known because A'hi is not measured in GRB 100901A). 
The two parameters are scaled until they reproduce the observed 
heavy-element depletion pattern best (minimum ;^f-). 

Results using measurements of 7 elements are shown in Fig- 
ure |4] The WH depletion pattern gives the smallest among the 
4 patterns, but the WD pattern is also a reasonable fit. The dust-to- 
metals ratio is between the MW value and 10% higher. Ni and Mg 
are the most uncertain measurements. The Mg ii column density is 
likely a lower limit because it suffers from very strong saturation. 
The Ni ii column density is not measured from the resonance lines 
directly, but follows from the modelling of the time evolution of 
the Ni n'''Fg/2 metastable level (see Section [3.2.1l and Table|5j. Ex- 
cluding these two elements from the fit does not change our results 
much. 

For the observed column of metals, and assuming a rate of 
visual extinction per column d ensity of metals like in the MW 
teohlin. Savage & Drake|[T978h . the expected optical extinction in 
GRB 100901 A would be Ay ~ 0.5. This is higher than the one de- 
rived from the spectral energy distribution (SED) of GRB 100901A 
(Ay = 0.21, Gomboc et al., i n prep). This incons istency is typically 
observed in GRB sight lines dSchadv et al.l201lh and indicates that 
the rate of extinction in the ISM of GRB hosts per column of metals 
is different from what is observed in the MW. 



3.1.3 Intervening absorbers 

Two intervening absorbers are detected: intl at zi = 1.3147 and 
int2 at zi = 1.3179. If these two systems are not physically as- 
sociated with each other, and the redshift difference is dominated 
by cosmological expansion, the velocity difference corresponds to 
a co-moving separation of 3.9 Mpc. If the absorbers belong to one 
system, the velocity difference is 414kms"', and could be due to 
motion of galaxies within a cluster. 

Because of the long time span that our observations cover 
(until long after the jet break, Gomboc et al., in prep), the data 
set is suited to look for variations in the strength of the lines 
from the intervening absorbers. This can be interesting in the con- 
text of the previously supposed discrepancy between the number 
of strong Mgii absorbers (IV,,(2796) > lA) in QSO and GRB 
sight lines; in the latter the re dshift number dens i ty was found 
to be about two times as hi gh OProchter et al. | 2006|^ Vergani et alj 
l2009l : ICucchiara et al.l2009l . However. ICucchiara et all ( 1201 2|) used 
a larger sample and did not reproduce this discrepancy. 

We do not detect variability above the 2cr level in any of 
the two intervening absorbers of GRB 10090IA, and we see no 
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Figure 4. Depletion pattern for the host of GRB 100901 A, compai ed with 
different components in the Milky Way jSavage & Sembachlll996ll . which 
are observed patterns that can be adjusted by changing the dust-to-metal 
ratio K. The column density of Ni ii is not directly measured but follows from 
the modelling of the fine-structure and metastable levels, see Section l3.2.1l 
We include a 3o" upper limit for Ti ii based on the absence of the doublet at 
a rest wavelength of 1910 A. 
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Figure 5. Curee-of-growth for the two intervening systems in the sightline 
to GRB 100901A. The lines used and the best-fit column densities are listed 
in Table|3] 
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systematic trend in the ensemble of lines per absorber, nor per ion 
species. In principle, the non-variation of the equivalent widths of 
the intervening absorber resonance lines gives a lower limit on the 
size of the absorbing clouds (i.e., the scale on which the absorbing 
gas is homogeneous). If the projected apparent size of the afterglow 
became larger than the projected size of the intervening absorber, 
the absorption lines from this system would become weaker since 
part of the light would reach the observer unabsorbed. We use the 
description of the appa rent size of the Blandford-M cKee spherical 
expansion described bv lGranot. Piran & Sari I(ll999l) . where we use 



6.3 X 10 erg dGorbovskov et al 



2012h . 



We assume the 

close circum-burst medium to have a constant density n = I cm"^^. 
At the time of the last epoch (7.027 days), we find the source to be 
9.6 X 10'* cm. This gives a lower limit on the size of each of the 
absorbing clouds of 0.03 pc. Although this is not a very meaningful 
lower size-limit for a gas cloud, we point out that in general it 



is diffic ult to put limits on absorber sizes, but s ee iD'Elia et al. 



2010al) for methods applied to GRB absorbers andlPetitiean et al 
20001): lEllison et 



jos apptiea to OKc absorbers andlhetitiean et al 
"ai] | |2004|) : iBalashev et al] l l20Ilh for examples 



of size estimates of QSO line-of-sight objects. 

Because the lines of the intervening systems do not vary 
in strength, we can average the equivalent widths measured in 
epochl and epoch2. The MISC-COG analysis can be applied, 
because we detect many clear lines from these systems. Figure |5] 
shows the best fitting COGs for both intervening absorbers. We 
carry out the same approach as for the host, by first constraining b 

for inti and 



with the Fe n transitions. We find b 



I2.6*\ j km s" 



b = 20.6^^ J km s ' for int2. The equivalent widths and the column 
densities of each ion species are listed in Table|3] 



3.2 Detection and variability of transitions arising from 
excited levels of Fe n and Ni n. 

At the redshift of the host galaxy we detect lines from fine-structure 
levels of Fen ^05/2, *£^3/2 and and metastable levels 

of Fen ('^Fg/2 and '^Dj^) and Nin (''^9/2). The equivalent width of 
many of these lines varies with time. This provides the opportunity 
to derive the distance between the burst location and the absorb- 
ing material, assuming that the population of these excited states is 
due to the UV-radiation of the afterglow (see Section l3.2.1) . With 
high-resolution spectra it would be possible to directly measure the 
column densities of the ions that are in the exci ted states, via Voigt 
profile fits to the fine -structure lines (see e.g., IP'Elia et alj|2007l ; 
IVreeswiik et al.ll2007h . In order to obtain column densities from 
the equivalent widths, we need to take an intermediate step with 
the COG, where we fix b to the value found from the resonance 
Fell lines (see Section [3.1. Il l: i.e., we assume that the ions in the 
excited states are in the same absorbing clouds as the ions in the 
ground state; for the implications of this assumption see the end of 
Section [3.2.1l Note that the value of b will only have a very small 
influence on the column densities of the excited ions, because the 
lines are all weak and lie mostly on the linear part of the COG. 

Tables IaTI I A2I I A3I and l A4l in the Appendix give an overview 
of all transitions from the fine-structure and metastable levels of 
Fe n and Ni 11 that would in principle be observable in the GMOS 
spectral range at this redshift. Due to the low spectral resolution, 
many of these lines are blended such that the individual variation 
of the components cannot be measured. The lines that are used in 
the modelling are shown in boldface in Tables |AT|to|A4] 

In order to obtain the column density Wievd of ions in a specific 



excited state in an epoch, the equivalent widths of the transitions 
arising from this excited level are placed on the COG by assign- 
ing an Wievci to this level. The best-fit A'lcvei is the value for which 
the ensemble of transitions fit the COG best (minimisation of the 
X')- Figures |6] and [7] show the COG of the Niii''F9/2 and the Feii 
''Dt/2 level per epoch (columns), together with the observed transi- 
tions from this level and gaussian fits to the line profiles. The best- 
fit A'levci and its errors are found with a Monte Carlo simulation, 
which we apply as follows. In every iteration, for each transition 
an equivalent width is randomly picked from a normal distribution 
with as mean the measured value and as sigma the measured er- 
ror. For this simulated ensemble of equivalent widths, the best-fit 
A'levci is determined by minimising Ihe x^- After 10000 iterations, 
the distribution of best-fit level's is fit with an asymmetric gaussian 
(i.e., a normal distribution with a diflterent cr on both sides), from 
which we derive the overall best A'levei and its lower and upper Itr 
error. Because the fine-structure lines lie mostly on the linear part 
of the COG, the effect of b is small, and the uncertainty of this value 
is not taken into account. If two or more transitions of a level are 
detected in an epoch (i.e., if > cr|Wi|), we obtain a value for 
A'lcvei, otherwise we obtain an upper limit (see Table|4j. 

3.2.1 UV-pumping models 

The column densities of the fine-structure and metastable states are 
found to vary in time (see Tableland Figure |8). The presence, 
and variation, of fine-structure lines in GRB afterglow spectra can 
generally well be explained by excitation due to th e UV flux of the 
GRB afterglow. IProchaska. Chen & BloomI ( l2006h suggest that this 
is the dominant excitation mechanism, especially when variability 
is measured. The early high-resolution spectroscopic data obtained 
for GRB 060418 jVreeswiik et al.l2007l) allowed systematic tests of 
other excitation mechanisms such as a background IR field and col- 
lisional excitation, but also strongly favours the UV-pumping sce- 
nario. For GRB 100901 A we assume that UV-pumping is the only 
relevant photo-excitation mechanism. 

We apply the photo- e xcitati on model introduced by 
IVreeswiiketaLl bOOTi 1201 iL |2013|) . For technical details we 
refer to the 2013 paper. The general idea is that the UV-flux of 
the afterglow temporarily excites ions in a cloud with thickness 
/ at a distance d from the location of the burst. As the afterglow 
brightness fades in time, the excited levels depopulate after a time 
which is determined by the Einstein coefficient A,,/ or, equivalently, 
the oscillator strength of the transitions considered. For a given set 
of input parameters, the column density of each excited level as a 
function of time, A'ievci(0. is predicted by the model. By comparing 
this temporal behaviour to the measured A'lcvei (0 for all levels 
simultaneously, we can optimise the parameters. 

The afterglow flux is included as an inteipolated series of R- 
band magnitude^ "'r(0 (see upper panel Figure |8), which corre- 
sponds to the rest-frame UV-flux (at ~ 2700 A) responsible for pop- 
ulating the levels. The monochromatic flux in the host-galaxy rest 
frame at the GRB -facing side of the absorbing cloud is computed 
as follows: 



l+Zh 



/lrcst(l + Zh) 



6415 A 



(3) 



The observed S-band magnit udes come from GCNs 



Andreev. Sergeev & Pozanenko* feOlOa'b'c'd'); lAndreev et al.l )201' 



Kuroda et ah i.2010a .b.): .Volnova et al. 1 2010) : assembled by Gomboc et al. 
(in prep.). 
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Table 3. The rest-frame equivalent widths for the lines from two intervening systems intl and int2 at respectively zi = 1.3147 and Z2 = 13 179. < W,{/A > 
is the weighted average of the rest frame equivalent width (in Angstrom) measured in GMOS epochl and epoch2. The COGs for the two systems are shown 
in Figure|5] The column densities per ion species have been calculated with the COG with the best fitting b; the en'ors on the column densities are the result of 
taking into account the range in b. Especially for lines on the flat part of the COG (e.g., Mg ii), the effect on N due to b is large. 







intl (71 
11111 


= 1.3147) 


int2 (z2 


= 1.3179) 






b = 12.6 




b = 20.6 


+5;^kms-' 


Line 


7,1 






<. VV j t\ > 


iog("ion'™ ' 


Aim 1854.7164 
Aim 1862.7895 


0.539 
0.268 


0.053 ± 0.046 
0.073 ± 0.047 


12.74^29 


0.229 ± 0.049 
0.171 ± 0.049 


13 42+"-27 
-0.17 


Fen 2249.8754 


2.190 X 10-3 


0.032 ± 0.036 




a 




Fen 2344.2129 


1.252 X 10-' 


0.317 ±0.030 


0.176 ±0.027 


13.57^°^ 


Fen 2374.4604 


3.297 X 10-- 


0.182 ±0.028 




0.034 ± 0.024 


Fen 2382.7641 


3.432 X 10-' 


0.365 ± 0.030 




0.325 ± 0.029 




Fen 2586.6495 


7.094 X 10-2 


0.296 ± 0.027 




0.136 ±0.023 




Fen 2600.1722 


2.422 X 10-' 


0.405 ± 0.028 




0.344 ± 0.027 




Mgn 2796.352 
Mgn 2803.531 


6.123 X 10-' 
3.054 X 10-' 


0.532 ± 0.030 
0.531 ± 0.030 


14.86+»;2^ 


0.845 ± 0.035 
0.792 ± 0.034 




Mgi 2852.9642 


1.83 


0.143 ±0.023 


12 17+°"' 


0.118 ±0.022 


12.01!JJ:?? 



"this line is blended with Ni n F9/2 /i2166.23 at the host redshift. 



Table 4. Lower level column densities log ^A'levciCO/'^rri for Fen and Nin excited states derived from fine-structure line equivalent widths with COG 
analysis. All errors and limits are Icr. 



Fe II Ni n 

time(d) ^Dgn ^Bia *D5/2 ^D^p_ *Dii2 ^Fg/^ 



epochl 0.0526 15.20!0;«8 13.77+«;05 13.58!0;« 13.30+«;0^ 12.92+«;li 12.54!";" 13.25!«;|« 

epoch2 1.0095 15.29!°;?^ 13.68l°;°5 13.46+°;»« 13.14+0;0» < 12.83 12.54+55;^] 13-V3+«;°^ 

xsh 2.7492 15.16+° [2 < 13.46 < 14.02 < 13.47 < 13.40 < 13.03 13.54+55°^ 

epoch3 7.0270 15.34;°-2l < 13.37 < 13.80 < 13.39 < 13.47 < 12.97 < 13.75 



in which Fo = 3.02 x IQ-^" erg s-' cm-^ Hz-' is the flux 
of Vega at the effective wavelength (6410 A) of the R-haxiA 
dFukugita. Shim asaku & Ichikawa" 1995) , As„ai = 0. 264 is the 
Galactic extinction ( Schlegel, Finkbeine r & Davisll998h . /l,cst is the 
wavelength amy of the relevant transitions at which the flux is re- 
quired to calculate the amount of excitation, is the spectral slope 
and Di = 9.968 x 10' pc is the luminosity distance. The effective 
Doppler parameter b is fixed to 22. 1 km s-' , consistent with how we 
converted equivalent widths to column densities (see Section [X2l l. 
Parameters that can be constrained from the light curve and/or SED 
fitting (Gomboc et al., in prep) are kept fixed: optical spectral slope 
jSy = 0.82, and optical extinction in the host galaxy Ay = 0.21 
assuming a Small Magellanic Cloud extinction profile (Gomboc et 
al., in prep). We performed models with Av = as well to study the 
effect. Ay is not included in Equation [3] because the model allows 
to specify where the optical extinction takes place. We place the ex- 
tinction in the absorbing cloud, and apply the necessary corrections 
to the flux. 

Table [5] shows the fit results of the model to the measured 
A'icvci(f) as listed in Table [4] with different requirements for the 
parameters. The model without extinction fits the data best with 
a small cloud at ~ 250 pc, but when we include optical extinc- 



tion Ay = 0.2, the data favours a configuration with a very large 
(> 1 kpc) absorbing cloud at a smaller distance. Fixing the cloud 
size / leads to absorber distances of rf ~ 165 - 275 pc. Figure [8] 
shows the fit to the data for the model with / < 500 pc (solid lines) 
and the model with / = 1 pc (dotted lines). Both models are decent 
fits: our data are not constraining enough to discriminate between 
the close-large and the distant-small cloud case. But in all cases, 
the influence of the GRB reaches distances of a few hundreds of 
parsecs. In Section|4T|we will discuss the implications of this es- 
timated absorber distance and compare it to values found in other 
GRBs. 

Though not detectable with the spectral resolution used, it 
could be possible that ground-state Fe 11 is present in more velocity 
components than excited Fe 11, i.e., that there is Fe n that is not asso- 
ciated with the excited gas. If there exists more Fe 11 further away, 
then this would decrease the amount of Fe 11 that is associated with 
the excited Feii, and hence this would move the absorber closer 
to the burst, i.e., we would infer a smaller GRB-absorber distance. 
Very roughly, if this fraction of Feii associated to excited Feii is 
only 50% (rather than the now assumed 100%), then it would de- 
crease the distance by a factor of V2. So even when assuming such 
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epoch2 
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Figure 6. Cui'ves-of-growth (COG) and spectram excerpts fortlie tliree detected lines from tfie Niii ^Fgi2 metastable level. Every column represents an epoch, 
with the observed time since the burst indicated on top. The upper plot shows the COG for b = 22.1 km s"' to which we fit the equivalent widths for the three 
transitions (red diamonds). The best-fit column density Wievci is indicated, together with the of the COG-fit. Below the COG, the three individual lines are 
displayed on a velocity scale, with the corresponding gaussian line profile fit. The rest frame equivalent width and oscillator strength are indicated. We use the 
same line label colouring as in Figure[T] 



Table 5. Photo-excitation modelling results for the environment of GRB 100901A with different settings. Column (1) extinction Ay in the host galaxy, (2) 
requirement for the distance d from the burst to the absorber, (3) requirement for the thickness / of the absorbing cloud, (4) fit result d, (5) fit result /, (6) 
pre-burst column densities of Fe ii and (7) Ni ii (ground state) and (8) the reduced chi-square of the fit. See also Figure[8] 

input output 

Ay d I d/pc l/pc log(A'Fen/cm"^) log(A'Niii/cm"2j xl 

0.0 free free 250 ± 75 Ot^*° 1^.25+55°^ l^.gi^jj jj^ 1.54 

0.2 minimal lOOpc free lOOt^' 1127 ± 569 15.26+55°^ 14.14!° |!J 0.73 

0.2 free fixed to 1 pc 275 ± 12 1 15-24+j5o^ 13.94+5J[j5 1.35 

0.2 free fixed to 100 pc 249 ± 13 100 15-24+°o6 13.95+°[55 1.33 

0.2 free maximal 500 pc 165 ± 73 500+ 15.22;[5 05 14.02+[' [!' 1.15 
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epoch 1 

0.0526 days after burst 



epoch2 

1 .0096 days after burst 



epochs 

7.0270 days after burst 




-500 500 
V (km s"^) 



-500 500 

V {km s^) 



-500 500 
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Figure 7. The same as Figure |6] only for the six detected lines from the 
Fen fine-stracture level. We omit the xsh epoch, because for this 

level none of the Unes was significantly detected due to the low S/N. 




D> 13.0 
o 



0.001 0.010 0.100 1.000 

1^1,3 (days since GRB trigger) 



10.000 



Figure 8. Top panel: observed S-band magnitudes from various observa- 
tories (see text), which are used as input for the photo-excitation model 
(see Section [3.2. It . Middle and bottom panel: column density as a function 
of observed time for Fe ii fine-structure and metastable levels, and for the 
Niii''Fc)/2 metastable level. The open circles and squares (arrows) indicate 
the column densities (limits) derived for each excited state (see Table |4). 
The solid lines are the predicted best fit result of the photo-excitation model 
with a maximum absorber thickness / < 500pc (xl = 115); the dotted line 



for a model with fixed absorber thickness / = 1 pc (xl 
Av = 0.21. 



1.35), both with 



a large fraction (50%) of the ground-state Fe n is further away does 
not change the infen'ed distance much. 

3.2.2 Atomic data 

In the full analysis from the COG to the photo-excitation model fit- 
ting, we need to adopt a set of values for vacuum wavelengths A and 
oscillator strengths /i. The results of the analysis described in this 
paper are obtained with the atomic data given in Tables lAniA2llA3l 
and IA4I This atomic param eter set is the same as the one used in 
CLOUDY dperland etal. Il998h. a widely -used photo-ionisation code, 
and c omes from IVemer et al.l (Il999li and references therein (see 
also iLedoux et al ]2009l) . The choice for this set is in principle 
arbitrary, but chosen such that our results can easily be compared 
with other works on GRB afterglow spectra. We tested the depen- 
dence of our results to this input by carrying out the full analysis 
with a different set of values for A and f,\. This alternative set is 
the same as used in 



Vr eeswiik et al. |2007h a nd uses a tomic data 

from lOuinet. Le Doumeuf & Zeippenl ( Il99a) . iMortonI ( l2003h and 
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1.190x10 1.195x10 1.200x10 1.205x10 1.210x10 

Wavelength (A) 

Figure 9. [Oiii]host emission lines, redshifted to ~ 1.2/j in the xsh-NIR 
spectram. The top panel shows the 2D spectram (0.5A/px). The bottom 
panel shows in black the normahsed ID spectram (0.5A/px); in blue the 
same spectram, but with a 10 pixel median filter to make the lines easier 
to see. The non-identified peaks in the ID spectram are due to residuals of 
skyline subtraction which are visible in the 2D spectram as vertical white 
strips. 



This leads to a 2% lower b, and differences up to 
0.06 dex in log N{t). The final resulting / and d of the cloud from the 
modelling in Section [3.2. ll differ by < 10%, sometimes with lower 
values fov xl- This latter fact is, however, not a reason to adopt the 
alternative set but rather to interpret the scatter in the results as an 
additional error. 



3.3 Emission lines 

Because GRBs originate in star forming galaxies, we have searched 
for emission lines at the host-galaxy redshift. There is tentative ev- 
idence for [O ni] /15007 (see Figure |9), of which we estimate the 
flux to be 4.2 ± 1.3 x 10"'^ erg s"' cm"^, error including systemat- 
ics due to flux calibration. [O iii] /14959 appears to be present, but is 
blended with a nearby skyline and cannot be measured. We do not 
detect Ha, H/? or [O n] due to low S/N and/or skyline blends. 



4 DISCUSSION 

We have shown that the GRB 100901 A afterglow UV radiation ex- 
cited Fe II and Ni ii ions in an absorbing cloud located a few hun- 
dred parsecs away in the host-galaxy ISM (Section|32}- The ionic 
column densities we derive for the part of this host galaxy that is 
along the line of sight (Section [3. It are very comparable to other 
GRB hosts wher e fine-structure lines are detected in the afterglow 
spectra (see e.g., PennraseetaL 20061: Prochaska. Chen & BloomI 



l2006l : IVreeswiik et al .120071: IP'Elia et al.l2009bh . Furthermore, the 
part of the galaxy we probe does not appear to be a special case 
in the context of dust depletion (Section [3. Lit . Modelling of the 
variation of the lines from the excited states allows us to estimate 
the distance to, and the size of, the absorbing cloud. This has been 
possible for about a dozen of bursts; in the following section we 
will discuss what these absorber distances teach us about the use of 
GRB afterglows as probes of galaxies. 



4.1 GRBs as probes of galaxies 

GRB afterglows are used as probes of their host galaxies, mak- 
ing use of the fact that the afterglow can be extremely bright for 
a short time. An advantage (for example over galaxies in quasar 
sight lines) is that afterglows select star-forming galaxies, under 
the generally accept ed assumption that GRBs are lin ked to (mas- 
sive) star formation (IWooslevlll99"3l: |PaczvnskillT998l) . In order to 
judge how useful they are as galaxy probes, it is important to un- 
derstand what the influence of the GRB and its afterglow is on the 
environment, and what part of the host galaxy we are probing. For 
about a dozen GRB afterglows (see Table |6), the distance between 
the region where the absorption lines originate and the site where 
the burst has gone off has been derived with spectroscopy. The typ- 
ical absorber distance is of the order of a few hundred pc, so the 
values we find for GRB 100901 A appear to be normal (Table |6j. 
These values are also in agreement with lower limits on the ab- 
sorber distance ( > 50 - 100 pc) revealed by the p resence of Mg i 
absorption lines (IProchaska. Chen & Bloomll2006h . 

Evidently, a GRB is able to affect the ISM out to distances of 
100 to 1000 pc, which can be much further out than only the star- 
forming region of the GRB progenitor. For example, the largest 
star-forming region in the Local Group, the Tarantula nebula, has a 
diameter of ~ 200 pc ( Lebouteiller et al. 20o3). Firstly, this means 
that the GRB and its afterglow could have considerably influenced 
the environment it probes. Secondly, the absorption lines that we 
observe are probably not from the immediate star-forming region of 
the progenitor of the burst. This is important, as many stellar evo- 
lution models that are proposed as long GRB progenitors critically 
rely on a low metallicity (in particular iron group element abun- 
dances, see e.g., IWooslev & Hegej|2006l) . Knowing the distance to 
the abundance-providing cloud gives more quantitative input into 
the reliability of the found abundances as tracers of the progenitor 
abundances. 



4.2 Can \o-w-R spectra be of use in the analysis of 
fine-structure line variation? 

Our analysis of the fine-structure line variability in GRB 100901 A 
differs from other studies because it is predominantly based on 
low-resolution spectra, while most other absorber distances are 
determined from spectra with a 25 - 40 times higher spectral 
resolution (see Table |6ll. Low resolution allowed us to moni- 
tor the spectrum over a much longer time-span. We note that 
GRB 100901 A remained particularly bright even after a week. The 
combination of our observations and the excitation modelling 
shows that we are detecting the peak of the Niii level pop- 
ulation at around 1 day after the burst in the observer's frame (see 
Figure [8}, which to date has not been observed. The downside is 
that low-resolution spectroscopy yields less accurate results, be- 
cause additional steps and assumptions are required to translate 
equivalent widths to column densities. Furthermore, it is difficult 
to disentangle different velocity components, limiting the degree 
to which we can map the internal dynamical structure of the host 
galaxy. 

When a (multi-epoch) low-resolution data set is available, we 
have shown that it is worthwhile to measure the variation and es- 
timate the absorber distance by fitting a photo-excitation model to 
it. Especially the lines from the Nin *Fg/2 metastable level have 
proven to be useful for these kind of long time-span spectroscopic 
data-sets. These lines are strong, relatively isolated from other com- 
mon lines in the spectrum, but still close to each other (with GMOS, 
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one would be able to catch all four of them if the burst is at 
0.8 < z < 1.9), and it takes longer for the levels to depopulate 
compared to the Fe ii fine-structure levels. In terms of recommen- 
dations for observing strategies, it would still be most useful to ob- 
tain a high-resolution spectrum as soon as possible after the burst, 
but late-time (after a few days depending on the brightness) spec- 
troscopic follow-up at low resolution can provide additional con- 
straints on the line variation. The information about the different 
velocity components and their Doppler parameters obtained from 
the earliest spectrum can be used to interpret the late time low- 
resolution data more accurately. 



5 SUMMARY AND CONCLUSIONS 

We have analysed the optical to near-IR spectra of the afterglow of 
GRB 100901A. Our data set consists of three low-resolution spec- 
tra obtained with Gemini-N/GMOS at 1 hour, 1 day and at seven 
days after the burst, and one medium-resolution VLT/X-shooter 
spectrum taken nearly 3 days post burst. At a redshift of 1.4084 
a wealth of metal resonance and fine-structure lines is detected, 
allowing us to estimate column densities of Fe, Si, Cr, Mn, Zn, 
Mg, Ca and Al with a multi-ion single component curve-of-growth 
(MISC-COG) analysis. The obtained metal column densities are 
similar to what has been found in many GRB sightlines. Ly-o- is 
not detected, so no metallicity can be determined. 

[Fe/Zn] = -1.13 ± 0.10 suggests the presence of dust in the 
host galaxy. A comparison of the o bserved dust depletion pat tern 
with sight lines in the Milky Way dSavage & Sembacl]|l99d) re- 
veals a slight preference for a warm-halo-like depletion pattern. 

The spectrum shows two intervening absorbers at redshifts 
1.3147 and 1.3179. Lines of Feii, Alii, Aim, Mgi and Mgii have 
been detected and the corresponding column densities are esti- 
mated with the MISC-COG technique. We do not see significant 
variation in the strength of any of these lines from the intervening 
absorbers. Under simple assumptions for the projected source size 
as a function of time, this lack of variability sets a lower limit of 
0.03 pc on the size (or homogeneity scale) of the absorbing clouds 
in these foreground galaxies. 

At the host-galaxy redshift we detect lines arising from fine- 
structure levels of the ground state of Fe n and metastable levels of 
Fen and Nin. The strengths of these lines are found to vary sig- 
nificantly in time. By grouping transitions that arise from the same 
lower level, we have used the COG from the Fe ii resonance lines 
to retrieve the temporal behaviour of the column density of each 
of these excited ion species. This information is fed into a photo- 
excitation model, that uses the rest-frame UV-flux of the burst as 
input. We assume that the excited levels are populated due to the 
UV-radiation produced by the afterglow. The model computes the 
distance between the burst location and the cloud of gas in which 
the lines originate, and the size of this cloud. For GRB 100901A, 
we estimate an absorber distance of a few hundred pc, which ap- 
pears to be typical compared with similar studies. 

This is the first time that a fine-structure line variability anal- 
ysis is predominantly based on low-resolution spectra, and as a 
consequence of that, also the longest time span over which such 
variation is detected and modelled. The decrease of the Niii"'F9/2 
metastable level population has not been detected before; our data, 
in combination with our modelling results, show that we cover the 
peak of the excitation of this level. We argue that applying the fine- 
structure variability model to low-resolution data can yield a sensi- 
ble estimate for the absorber distance, but additional intermediate 
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steps (such as the COG) are necessary, which may cause additional 
(systematic) errors. 
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Table 6. Bursts for which an absorber distance has been derived from (the variability of) lines from excited and/or ionised states in the afterglow spectrum. 
Column (1) Burst ID, column (2) number of epochs, column (3) approximate time (observer's frame) between the different epochs, column (4) telescope and 
spectrograph with which the data have been taken, columns (5) and (6) the ions and levels/lines for which variability is measured and/or modelled, column (7) 
redshift of the burst, column (8) neutral hydrogen column density columns (9) and (10) distance between the burst and the region where the lines originate, 
and size of this region following from modelling of the variability, column (11) reference. 



Burst 



# 



instrument 



ion levels/lines 



logWffl 



J/pc 



//pc 



ref. 



020813 2 

050730 2 

051111 1 

060206 2" 



16h 
Ih 



KECK/LRIS 

VLT/UVES 

VLT/UVES 



KECK/HIRES 
WHT/ISIS 



060418 6 5 -30min VLTAJVES 



080310 2-4 10-20min VLTAJVES 



0803 19B 3 40-60min VLTAJVES 



080330 1 



VLTAJVES 



081008 4" 7-25min VLTAJVES 

VLT/FORS2 
090926 4" VLT/X-shooter 



100219A 1 VLT/X-shooter 
100901A 4 1 h - 1 wk Gemini-N/GMOS 

VLT/X-shooter 



Fell *D7/2 /12396 



1.25 



Veil ^Dii2,^Dii2°Dil2°Dii2, 3.97 22.10 



*f9/2, Ft/2,F5/2,F3/2, 



"01,2, *Di,2 
Fen *D7/2, *I>5/2,*£'3/2,''£'i/2, 



3/2 



Sill 

Sin 2p°/2, 
Oi 3P°,'P° 



1.55 



4.05 20.85 ±0.1 



50 - 100 



124 ± 20 



a few times 10 



147+ 



Cii 



3/2 



1.49 



10^ 



480 ± 56 



Fell "0^2,05/2,03/2,01/2, 

^Fg/2, *Dt/2, 
Nin '^Ftj/2 

Fell *Z)7/2,'Z)5/2,*£>3/2.'^Oi/2, 2.43 18.70±0.1 200-400 
'*Fc)/2, ''D7/2, 



0- 100 



Fern 'Z), '53, 
Sin 

3/2 



Cii 



3/2 



Fen '■Z)7/2, ''Z)5/2, O3/2, O1/2, 

''Fq/2, *Ft/2, 
*Oi/2,*Oii2, 
Fen ^Oi/2,^Oi/2,^03/2,^Oi/2, 

*Fq/2, 

Sin '-Py, 
Nin '^Fg/2 

Fen ^Ot/2,^Os/2,^03/2. 

*Fq/2, *Ol/2, 

Sin 

Nin ''F9/2 

Fe n ^Ot/2, '*Fg/2, 

Sin 

3 DO 3, 
Cn 

Nin *Fg/2 
Sin ~Py2 

Fen ^Ot/2,^Os/2,^03/2,^Oi/2, 

*Oj/2, 
Nin '*Fi,/2 



0.94 



1.51 



560 - 1700 



1.97 21.11 ±0.10 52 ± 6 and 200+f° 



2.11 21.60 ±0.07 677 ± 42 and 5 X 10^^ 



4.67 
1.41 



21.14±0.15 



300 and lO' 
a few times 10^ 



1 - 1000 



2,3,4 



5,6 



7, 8 



9,3 



10, 11 



12, 3 



13, 3 



14 



15, 3 



16 
17 



"More than one epoch is obtained, but the distance measurement did not follow from modelling the vaiiation. 

References: f 1) Dessauges-Zavadskv et al. (2006) : (2) Ledoux et al. (2009) : (3) Vreeswiik et al. (2011) : (4) D'Elia et al. (2007) : (5) Penprase et al. (2006) : 
(6) Prochaska. Chen & Bloom (2006); (7) Fvnbo et al. (2006); (8) Thone et al. (2008); (9) Vreeswiik et al. (2007) : (10) De Cia et al. ( 2012) : 
(11) Vreeswijk et al. (2013); (12) D'Elia et al. (2009a) : (13) D'Elia et al. (2009b) ; (14) D'Elia et al. (20 iT) ; (15) D'Elia et al. (2Q10b) : 
(16) Thone et al. (2012) : (17) this work. 
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Table Al. Lines from the ground state of Fen (electron configurations 3rf*('D)4.s- 3d^(^D)4p) and its fine-structure levels, that could in principle be detected 
in the spectrum of GRB 100901 A, because the transition would lie within the GMOS spectral range and the oscillator strength f,\ is not smaller than about 
0.005 (except for the ^Dgji resonance lines). The lines of which the wavelengths are in boldface have been used in the modelling reported in Section |X21 Wa 
are measured equivalent widths per epoch, converted to the rest frame. 



lower level 


'Ivac (A) 


A 


logs/ 




Jk 


terms 




W,i(epochl) 


W,i(epoch2) 


H',i(xsh) 


M',,(epoch3) 




1608.4509 


0.05399 


-0.268 


9/2 


-7/2 


a^D- 




possibly present, but blended with Al n /il670 at z 


= 1.318. 


cm-' 


1611.2004 


0.00136 


-1.867 


9/2 


-7/2 


a^D- 




not detected. 










2249.8754 


0.00219 


-1.660 


9/2 


-7/2 


a^D- 


z'D" 


0.129 ±0.034 


0.151 ±0.040 


0.143 ±0.044 


0.093 ±0.138 




2260.7793 


0.00262 


-1.582 


9/2 


-9/2 


cfiD- 


,4po 


0.152 ±0.034 


0.168 ±0.040 


0.074 ± 0.060 


0.256 ±0.125 




2344.2129 


0.12523 


0.098 


9/2 


-7/2 


a^D~ 




0.748 ± 0.047" 


0.798 ± 0.054" 


0.704 ±0.189" 


0.557 ±0.134" 




2374.4604 


0.03297 


-0.482 


9/2 


-9/2 


a*D- 




0.590 ± 0.045 




b 


b 




2382.7641 


0.34321 


0.536 


9/2 


- 11/2 


a«Z)- 


Z^F° 


0.683 ± 0.045 




b.c 


b 




2586.6495 


0.07094 


-0.149 


9/2 


-7/2 


a^D- 


Z^D" 


0.705 ± 0.043 


0.727 ± 0.048 


0.696 ±0.145 


0.629 ±0.107 




2600.1722 


0.24226 


0.384 


9/2 


- 11/2 


a^D- 


z'D" 


0.780 ± 0.045 


b 


b.c 


b 



"Blended with lines from 'Z)i/2 and '^F-ij2, total W,{ is reported. 

*Wj of this resonance line is kept fixed to the value obtained at the first epoch to see the possible variation of the (half-) blended fine-structure lines. 
'^The line (blend) is clearly saturated, and cannot be fit with gaussian (components) in the X-shooter spectrum. 



«D7/2 2333.5147 


0.07776 


-0.206 


7/2 


-5/2 


a^D- 


,6po 


0.168 ±0.034 


0.137 ±0.038 


0.051 ±0.189 


-0.134 ±0.126 


384.790 cm-' 2365.5508 


0.05285 


-0.374 


7/2 


-7/2 


a^D~ 


z'P° 


0.185 ±0.036 


0.130 ±0.039 


0.008 ± 0.222 


0.056 ±0.120 


2383.7873 


0.00557 


-1.351 


7/2 


-5/2 


a^D~ 


z'F" 


0.087 ±0.031 


0.071 ±0.034 


C 


0.054 ±0.107 


2389.3569 


0.08987 


-0.143 


7/2 


-7/2 


a^D~ 


,6po 


0.187 ±0.033 


0.216 ± 0.040 


-0.085 ± 0.303 


-0.056 ±0.112 


2396.3551 


0.30779 


0.391 


7/2 


-9/2 


cfiD- 


z^F" 


0.383 ± 0.040'' 


0.348 ± 0.043'' 


0.170 ±0.275'' 


0.314 ±0.114'' 


2599.1457 


0.10863 


-0.061 


7/2 


-5/2 


a^D- 




0.190 ±0.031 


0.096 ± 0.030 


C 


-0.019 ± 0.096 


2612.6536 


0.12485 


-0.001 


7/2 


-7/2 


a^D- 


z'D" 


0.297 ± 0.035 


0.241 ±0.037 


0.012 ±0.106 


0.049 ± 0.096 


2626.4503 


0.04266 


-0.467 


7/2 


-9/2 


a^D- 


z^D" 


0.144 ±0.030 


0.120 ±0.031 


0.095 ±0.118 


0.012 ± 0.092 



''Blended with *i)5/2/12396. 1487, but we assume that the contribution from this line to the total W,t is negligible. 



«D5/2 2328.1100 


0.03760 


-0.647 


5/2- 


3/2 


a^D- 


z'P° 


0.052 ±0.031 


0.031 ±0.036 


0.064 ± 0.235 


0.082 ±0.132 


667.683 cm-' 2349.0215 


0.08935 


-0.271 


5/2- 


5/2 


a^D~ 




0.204 ± 0.035'' 


0.196 ±0.038'' 


0.081 ±0.1 54'' 


-0.139 ±0.114'' 


2381.4877 


0.03798 


-0.642 


5/2- 


7/2 


cfiD- 


,6po 


0.022 ± 0.027 


-0.009 ± 0.033 




-0.004 ±0.111 


2396.1487 


0.01636 


-1.008 


5/2- 


3/2 


d'D- 


z'F" 


blend with 


7/2/12396.3551, which dominates. 




2399.9718 


0.12522 


-0.124 


5/2- 


5/2 


a^D- 


Z^F° 


0.169 ±0.032 


0.163 ±0.037 


0.071 ±0.255 


0.221 ±0.124 


2405.6173 


0.25103 


0.178 


5/2- 


7/2 


d'D- 


Z^F° 


0.254 ± 0.036^ 


0.203 ± 0.039'' 


0.005 ± 0.205 ' 


-0.076 ±0.102' 


2607.8658 


0.12304 


-0.132 


5/2- 


3/2 


a^D- 


z'D" 


0.308 ± 0.040* 


0.242 ± 0.043 


0.083 ±0.108 


0.007 ±0.131 


2618.3984 


0.04800 


-0.541 


5/2- 


5/2 


ifiD- 


,6po 


0.112 ±0.029 


0.057 ±0.031 


0.127 ±0.109 


-0.038 ±0.114 


2632.1077 


0.08240 


-0.306 


5/2- 


7/2 




z^'D" 


0.252 ± 0.033* 


0.169 ±0.0358 


-0.015 ±0.114* 


0.086 ± 0.088* 



''Blended with *F9/2/12348.8346, this is the total Wa- 

^Blended with 'Z)3/2/12405.1626, this is the total Wa- 

*Blended with Mnn /12606, but this line's Wa is fixed after the first epoch. 

''Blended with '03/2^2631.8314, fills is the total Wa. 





2338.7238 


0.09840 


-0.405 


3/2- 


3/2 


a^D- 


-6po 


0.090 ±0.032 0.057 ±0.036 0.127 ±0.174 


0.062 ±0.128 


862.613 cm-' 


2359.8270 


0.06788 


-0.566 


3/2- 


5/2 


cfiD- 


,6po 


part of a strongly blended complex, with a contribution of an unidentified line. 




2405.1626 


0.02784 


-0.953 


3/2- 


1/2 




,(.po 


blend with ''05/2/12405.6173, which dominates. 






2407.3934 


0.15988 


-0.194 


3/2- 


3/2 


a^D- 


,6po 


0. 122 ± 0.032 0.078 ± 0.033 0.014 ± 0. 199 


0.014 ±0.112 




2411.2522 


0.21967 


-0.056 


3/2- 


5/2 


a^D- 


z'F" 


0.227 ±0.035' 0.115 ±0.037' 0.181 ±0.170' 


0.031 ±0.089' 




2614.6047 


0.11018 


-0.356 


3/2- 


1/2 


a^D- 


z'D" 


0.119 ±0.029 0.091 ±0.032 0.007 ±0.105 


0.055 ± 0.095 




2621.1906 


0.00364 


-1.837 


3/2- 


3/2 


a^D- 


z'D" 


0.047 ±0.028 -0.016 ±0.032 0.035 ±0.114 


0.022 ± 0.084 




2631.8314 


0.12430 


-0.303 


3/2- 


5/2 


J'D- 


z^D" 


blended with ''D5/2/12632.1077, see there for the total Wa. 



'Blended with ''Di/2/i2411.8009, this is the total Wa. 



^Dii2 2344.9996 


0.16819 


-0.473 


1/2- 


3/2 


a^D- 


,6po 


blended with 


'Z)9/2/12344.214 and ''F7/2/12344.6789, see there for the total Wa. 


977.053 cm-' 2411.8009 


0.22500 


-0.347 


1/2- 


1/2 


a^D- 


j6po 


blended with 


*Z)3/2/1241 1.2522, see there for the total Wa. 


2414.0437 


0.19047 


-0.419 


1/2- 


3/2 


a^D- 


z'F° 


0.060 ± 0.030 


0.028 ±0.034 -0.081 ±0.159 0.073 ±0.125 


2622.4513 


0.05454 


-0.962 


1/2- 


1/2 


a^D- 


z'D" 


0.036 ± 0.025 


0.024±0.029 0.114±0.109 0.116±0.114 


2629.0769 


0.17451 


-0.457 


1/2- 


3/2 


a^D- 


z^D" 


0.079 ± 0.027 


0.033 ±0.031 0.008 ±0.172 0.119 ±0.109 
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Table A2. The same as Table lAll only for the hnes from the metastable states of Fen with electron configuration 3p*3rf'. 



luwci lev CI 


'^vac 


f. 

J A 






J, 

•Jk 








"^-9/2 


1610.9234 


0.00914 


-1.039 


9/2 


-9/2 


a-*F- 




not detected, hardly in range, noisy region. 


1872.567 cm"' 


1637.3998 


0.01418 


-0.848 


9/2 


- 7/2 


a'^F- 


x*D° 


not detected, noisy region with many possible lines. 




1 709 OAS'? 


U.UJ J lO 


-u.z 


Q/9 

y/z 


1 1 /9 

— 1 1 /z 


CI r — 


z ^ 


possibly detected, but shifted and blended with Ni n/il703, noisy region. 




2332.0235 


0.01787 


-0.748 


9/2 


-7/2 


a*F- 


z*F° 


0.052 ±0.029 0.085 ±0.037 0.109 ±0.199 0.103 ±0.123 




2348.8346 


0.03653 


-0.437 


9/2 


-7/2 


a^F- 


z'D" 


blended with ''Z)5/2/12349.0215, see there for total Wx. 




2360.7213 


0.02289 


-0.640 


9/2 


-9/2 


a'^F- 


z*F° 


part of a strongly blended complex, with a contribution of an unidentified line. 




1625.5231 


0.02481 


-0.702 


7/2 


-9/2 


cdF- 


/G° 


not detected, noisy region. 


2430.097 cm-' 


2344.6790 


0.01669 


-0.874 


7/2 


-5/2 


a^F- 


z*F° 


blended with ''1)9/2^2344.214 and 'Z)|/2i2344.9996, see there for the total Wa- 




2361.0161 


0.03053 


-0.612 


7/2 


-5/2 


a^F- 


z'D" 


part of a strongly blended complex, with a contribution of an unidentified line. 




2380.0023 


0.01732 


-0.858 


7/2 


-7/2 


a^F- 


z'D" 


not detected. 




1633.9093 


0.02620 


-0.803 


5/2 


-7/2 


a^F- 


v*G= 


not detected, noisy region with many possible lines. 


2837.950 cm-' 


2369.3195 


0.02918 


-0.757 


5/2 


-3/2 


a^F- 


z'D° 


not detected. 




2383.9718 


0.02215 


-0.876 


5/2 


-5/2 


a^F- 


z'D" 


blended with many lines including the saturated *Dg/2/l2382.7641. 




2375.9188 


0.03512 


-0.852 


3/2 


-1/2 


a'F- 


v^D° 


0.061 ±0.034^ 0.066 ±0.034^ -0.147 ± 0.202^ -0.079 ± 0.100^ 


3117.461 cm-' 


2385.1148 


0.02311 


-1.034 


3/2 


-3/2 


a*F- 




not detected, close to complex blend. 



^Half-blended with *Z)9/2/12374.4604, which's W,, is fixed from the first epoch. 



Table A3. The same as Table lAll only for the lines from the metastable states of Feii with electron configuration 3d'(^D)4«. 



lower level 


-Ivac (A) 








Jk 


terms 




H',i(epochl) Wi(epoch2) W,i(xsh) H',,(epoch3) 




1635.4005 


0.06857 


-0.261 


112 


-5/2 


a'^D- 




heavily blended, noisy region. 


7955.299 cm"' 


2563.3039 


0.13520 


0.034 


112 


-5/2 


a'^D- 




0.037 ±0.027 0.064 ±0.032 -0.052 ± 0.248 0.047 ±0.110 




2715.2173 


0.04543 


-0.440 


7/2 


-5/2 


a*D- 


z^D" 


not detected. 




2740.3577 


0.25107 


0.303 


7/2 


-7/2 


a'^D- 


z'D" 


0.056 ±0.026 0.027 ±0.029 0.018 ±0.099 0.114 ±0.094 




2756.5512 


0.30759 


0.391 


7/2 


-9/2 


a'^D- 


z'F° 


0.050 ±0.025 0.056 ±0.030 0.097 ±0.114 -0.045 ±0.104 


"05/2 


1641.7631 


0.04742 


-0.546 


5/2 


- 3/2 


a*D- 


x*P° 


blended, noisy region. 


8391.938 cm-' 


1647.1625 


0.02026 


-0.915 


5/2 


-5/2 


a'^D- 


x^P" 


blended, noisy region. 




2564.2441 


0.09332 


-0.252 


5/2 


-3/2 




z'P° 


not detected. 




2592.3179 


0.05219 


-0.504 


5/2 


-5/2 


a*D~ 


z*P° 


not detected. 




2728.3467 


0.06890 


-0.384 


5/2 


-3/2 


a*D~ 


z*D° 


not detected. 




2747.7942 


0.18791 


0.052 


5/2 


-5/2 


a*D~ 




part of strongly blended complex. 




2750.1347 


0.32808 


0.294 


5/2 


-5/2 


a*D- 


z*F° 


blended with "03/2^2749.9938 and ''Di/2/i2750.2989. 


"03/2 


2567.6820 


0.05684 


-0.643 


3/2 


-1/2 




,4po 


not detected. 


8680.454 cm-' 


2583.3564 


0.08345 


-0.477 


3/2 


-3/2 


a*D~ 


j4po 


not detected. 




2731.5431 


0.02964 


-0.926 


3/2 


-3/2 


a*D- 


z*F° 


not detected. 




2737.7761 


0.06574 


-0.580 


3/2 


-1/2 


a*D- 


z'D" 


not detected. 




2747.2958 


0.35475 


0.152 


3/2 


-5/2 


a'^D- 




part of strongly blended complex. 




2749.9938 


0.13492 


-0.268 


3/2 


-3/2 


a'^D- 


z'D° 


blended with '♦05/2/12750.1344 and ''Z)i/2i2750.2988. 


"•01/2 


2578.6940 


0.12362 


-0.607 


1/2 


-1/2 


a'^D- 


z'P° 


not detected. 


8846.768 cm-' 


2744.0089 


0.44026 


-0.055 


1/2 


-3/2 


a'^D- 


z'F° 


possibly detected, but half blended Mgi /i2853 at z = 1.315. 




2750.2989 


0.12021 


-0.619 


1/2 


-1/2 


d^D- 


z'D" 


blended with ''05/2/12750.1347 and ''03/2/12749.9938. 



Table A4. The same as Table lAll only for the lines from the metastable states of Ni ii with electron configurations ^p^icfii^ F)As. 



lower level 


/Ivac (A) 


A 


log?/ 


/, - 


Jk 


terms 




W.i (epoch 1) 


W,i(epoch2) 


W,((xsh) 


M',i(epoch3) 


"^9/2 


2166.2300 


0.09990 


0.000 


9/2 


-9/2 


4f _2 


F° 


0.070 ± 0.035 


0.191 ±0.045 


0.197 ±0.046 


0.026 ±0.145 


8393.90 cm-' 


2217.1676 


0.17865 


0.252 


9/2 


- 11/2 


4f _4 


G° 


0.1 10 ±0.035 


0.222 ± 0.042 


0.150 ±0.043 


0.227 ±0.143 




2223.6419 


0.04789 


-0.320 


9/2 


-9/2 


4f _4 


G° 


0.019 ±0.031 


0.116 ±0.039 


0.042 ± 0.046 


-0.060 ±0.136 




2316.7481 


0.11073 


0.044 


9/2 


-7/2 


4f _4 


D° 


blended with 


a skyline. 
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